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The Extended X-ray Halo of the Crab-like SNR G21.5-0.9 
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Abstract. Recent XMM-Newton observations reveal an extended (r ~ 150") low-surface brightness X-ray halo 
in the supernova remnant G21.5-0.9. The near circular symmetry, the lack of any limb brightening and the non- 
thermal spectral form, all favour an interpretation of this outer halo as an extension of the central synchrotron 
nebula rather than as a shell formed by the supernova blast wave and ejecta. The X-ray spectrum of the nebula 
exhibits a marked spectral softening with radius, with the power-law spectral index varying from V = 1.63 ± 0.04 
in the core to T = 2.45 ±0.06 at the edge of the halo. Similar spectral trends are seen in other Crab-like remnants 
and reflect the impact of the synchrotron radiation losses on very high energy electrons as they diffuse out from 
the inner nebula. A preliminary timing analysis provides no evidence for any pulsed X-ray emission from the core 
of G21.5-0.9. 
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1. Introduction 

Of the ~ 225 known supernova remnants (SNRs) in our 
galaxy no more than 5% are classified as Crab-like rem- 
nants (Green 2000). Prominent members of this class of 
SNR, also known as plerions, include the Crab Nebula, 
CTB 87, and 3C58. The characteristic centre-filled radio 
and X-ray morphologies of the Crab-like SNRs is thought 
to be due to the presence of an active pulsar which powers 
a bright synchrotron nebula. In the case of the Crab neb- 
ula, 33 ms pulsations testify to the presence of the central 
spinning neutron star, although in many Crab-like systems 
there is no direct evidence for pulsed emission. Typically 
a Crab-like SNR exhibits a flat power-law spectrum in the 
radio regime which eventually steepens at shorter wave- 
lengths to join smoothly to a similarly featureless power- 
law X-ray continuum. This non-thermal form for the radio 
to X-ray spectrum is a consequence of the continuous in- 
jection by the pulsar of high energy electrons, which suffer 
radiation and adiabatic losses as they diffuse through the 
nebula; however the details of the process are undoubt- 



edly complex (e.g. Reynolds & Chanan 1984). A feature 
of the Crab-like remnants, which distinguishes them from 
both shell-like and composite-type SNRs (the latter typi- 
cally exhibiting a filled-centre X-ray morphology within a 
radio shell) , is the lack of any evidence for an outer shell 
structure marking the progress of the blastwave from the 
original supernova explosion. 

G21.5-0.9 is a SNR with many of the characteristics 



of a Crab-like remnant (e.g. Becker & Szymkowiak 1981 
Asaoka & Koyama 1990| ) . To date there has been no detec 



tion of pulsed emission in either the radio or X-ray bands 



Send offprint requests 
rsw@star .le.ac.uk 



to: 



R.S. Warwick, e-mail 



(e.g. Biggs & Lyne 1996; Slane et al. 2000). However, re- 
cent X-ray observations made by Chandra have pinpointed 
the probable location of the pulsar, namely a very com- 
pact central core on a scale of ~ 2" embedded within a 
more extended (~ 30" radius) synchrotron nebula (Slane 
et al. |2000 ). A faint extended halo was also detected which 
might correspond to the outer "shell" formed from the 
expanding ejecta and the pa ssage of the supernova-driven 
blastwave (Slane et al. 2000| ). In this interpretation G21.5- 
0.9 would be best described as a composite remnant. Here 
we report further recent X-ray observations of G21.5-0.9 
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Fig. 1. Left panel: The raw XMM-Newton EPIC PN image of G21.5-0.9. The pixel size is 4" and the image is 8' on 
a side. A logarithmic intensity scaling has been applied. The vertical and horizontal white bands correspond to the 
gaps between adjacent CCD chips of the PN camera. Right panel: EPIC MOS image of G21.5-0.9. The pixel size is 1" 
and the image is 8' on a side. The image has been spatial filtered with a Gaussian smoothing mask with width a = 2 
pixel. Again logarithmic intensity scaling has been applied to emphasise low surface-brightness features. 



made by XMM-Newton which help to clarify the nature 
of the extended X-ray halo in this object. 



2. XMM-Newton observations 

G21.5-0.9 was observed as a calibration target during or- 
bits 60-62 and 64-65 of the XMM-Newton mission. Here we 
focus on measurements from orbit 60 (date: 2000/04/07; 
start time: 12.35 UT; end time 22.57 UT) in which the 
SNR was observed on- axis b y the EPIC MOS a nd PN 
cameras (Turner et al. |200l| ; Striider et al. |200l| ) giving 
a total accumulated exposure time of ~ 30 ks. For these 



observations all three EPIC cameras were operated in the 
standard full-frame mode with the medium filter selected. 
Attitude information files are not yet available but all the 
indications are that a stable pointing was maintained for 
the duration of the observation^. 

The recorded events were screened with the XMM 
Science Analysis Software (SAS) to remove known hot 
pixels and other bad data and pre-processed using the lat- 
est CCD gain values. X-ray events corresponding to pat- 
terns 0-12 for the two MOS cameras (similar to grades 
0-4 in ASCA) were used, whereas for the PN only pattern 
events (single pixel events) were accepted. Investigation 
of the full-field count-rate revealed a number of flaring 
events in the non-cosmic component of the background 



1 Our nominal position for the centre of G21.5-0.9 is (RA, 
Dec) = 18 hr 33 min 33.8 s, -10° 34' 6" (J2000); (1, b) = 
21.50°, -0.88°. 



during the observation. As a final step in the data screen- 
ing we identified a period when the background was essen- 
tially "quiescent" , with full-field background rates in the 
MOS and PN cameras of ~ 2.5 and ~ 6 count/s. Selecting 
events only in this low-background period lead to effective 
exposure times of 13800 s and 11800 s for the MOS and 
PN cameras respectively. 



3. Spatial Analysis 

Given the observed soft X-ray cut-off in the spectrum of 
G21.5-0.9 (see §4), we utilize the 1-10 keV energy band in 
the spatial analysis. The raw image of G21.5-0.9 from the 
PN dataset binned into 4" pixels is shown in Fig. [l] {left 
panel) . The bright synchrotron nebula at the centre of this 
SNR is the most prominent feature of the image, although 
the central core is clearly surrounded by a more extended 
low surface-brightness emission region. In addition a point 
source is evident ~ 2' to the south-west of the centre of 
the remnant. 

To avoid the complications due to the chip boundaries 
in the PN data we have based our spatial analysis on the 
data from the MOS cameras (since the 10.9' x 10.9' field 
of view of the central chip in the MOS cameras comfort- 
ably encompasses the full extent of the SNR). Broad band 
images were constructed for both MOS 1 and 2 cameras 
utilizing 1" pixels. The MOS 1 and 2 images were then 
aligned (by eye) using the nebula core and the south- 
western point source as fiducial marks and co-added. Fig. 
[j] {right panel) shows the resultant image after applying 
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Fig. 2. The radial profile of the X-ray surface brightness 
of G21.5-0.9 (clotted curve). The counts in the broad-band 
MOS image were accumulated in annuli of 10" width out 
to a maximum radius of 200". The mean background level 
in the annular region from r — 200" to 240" was then sub- 
tracted. Note that within a radius of 150" the statistical 
errors are smaller than the plotted points. The presence of 
the south-western point source has a marginal impact at 
r = 115". By way of comparison, the point spread func- 
tion of mirror module 3 (MOS 1) at 1.5 keV i s also shown 
(dashed-dotted line) (Aschenbach et al. |2000|) . 



a Gaussian smoothing filter with a width a — 2 pixels 
(FWHM of 4.7"). At a qualitative level there is clearly 
excellent agreement between the images from the PN and 
MOS cameras. The bright centre of the synchrotron neb- 
ula (i.e. the region within r w 30" of the centre) shows 
deviations from circular symmetry, in particular there is 
an indentation in the north-west quadrant. At a slightly 
lower level of surface brightness a curious "spur" is ev- 
ident which appears initially to track northwards away 
from the core region but then at r ~ 100" sweeps round 
in an arc into the north-eastern quadrant of the rem- 
nant. These spatial features match rather well to simi- 
lar structure seen in the higher spatial resolution images 
from Chandra (Slane et al. 2000). One of the strengths 
of the EPIC cameras on XMM-Newton is, of course, the 
sensitivity afforded to extended, relatively low-surface X- 
ray emission. This is evident in the images of Fig. |l| in 
which a low surface brightness "halo" is well delineated 
out to r w 150". This halo exhibits a near perfect cir- 
cular symmetry apart from the hint of a "linear termina- 
tion" along its northern rim (note that out-of-time events, 
namely events recorded during the period when charge is 
transferred to the readout node, partly mask this feature 
in the PN image). 

In order to investigate more quantitatively the spa- 
tial extent of the low surface brightness halo, we have 
calculated the azimuthally-averaged radial profile of the 
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Fig. 3. The softness ratio image of G21.5-0.9 derived from 
the MOS 1 + MOS 2 datasets. Softness ratio is defined 
here as SR = (S-H)/(S+H), where S and H represent the 
1-3 keV and 3-10 keV band images respectively. The S 
and H band images were initially heavily smoothed with 
a tophat filter of width 30 pixels (1 pixel = 1"). Red 
to blue colours reflected a decreasing softness ratio. The 
black contours are based on the combined S+H smoothed 
image. 

X-ray emission (Fig. |2|). The surface brightness away from 
the centre of the synchrotron nebula drops rapidly out to 
r « 50". A plateau region is then reached corresponding 
to the low surface brightness halo. A sharp further decline 
in the surface brightness then sets in at r « 130", with 
the signal eventually falling below the detection threshold 
at r w 180". In this paper we define the remnant's outer 
radius to be 150" at which point the surface brightness is 
roughly one-fifth of the plateau value. 



4. Spectral Analysis 

We first examined the spectral properties of G21.5-0.9 by 
plotting the softness ratio image shown in Fig. ^, which 
shows that the core of the nebula is significantly harder 
than the low-surface brightness halo. 

We have investigated the spectral characteristics of the 
core of G21.5-0.9 using data from both the PN and MOS 
cameras. The X-ray spectrum of the core of the nebula 
was obtained using a circular extraction region of radius 
48" centred on the position of peak surface brightness. 
Background spectra were extracted from representative 
regions greater than 200" away from the SNR. The source 
spectra were binned to a minimum of 20 counts per spec- 
tral channel, in order to apply \ 2 minimisation techniques. 
The source and corresponding background spectra were 



4 



Warwick et al.: XMM-Newton Observations of G21.5-0.9 




Fig. 4. Upper panel: The observed PN (blue), MOS 1 (red) 
and MOS 2 (green) count rate spectra for the core re- 
gion of G21.5-0.9. In each case the corresponding best fit- 
ting spectral model is represented by the solid histogram. 
Lower panel: The ratio of the predicted to the observed 
count rate spectra. 

Table 1. Best-fitting spectral parameters for the SNR 
core. 
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o qo+0.05 
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1 8R+ - 04 


l.OD_ 07 
1 rc+0.09 
1 -°°-0.09 

1 61 +009 

1 -°- L -0.09 


1278/1203 
384/425 
448/423 



a Column density in units of 10 
6 Photon spectral index 

c Normalisation in units of 10 -2 photon cm -2 s _1 keV~ 



then analysed using the xspec vll.01 and the most re- 
cent response matrices available from the EPIC team. The 
spectral modelling assumed a power-law continuum (pho- 
ton spectral index, T, and normalisation, A) subject to 
soft X-ray absorption in cool gas along the line of sight 
(with an equivalent hydrogen column density, -/Vh)- 

Fig. ^ shows the observed count rate spectra for the 
PN, MOS 1 and MOS 2 detectors together with the corre- 
sponding best-fitting model spectra. The spectral response 
of the two MOS cameras is so similar that the MOS 1 and 
MOS 2 count rate spectra almost totally overlap. The fig- 
ure illustrates that the observed count rate in the PN cam- 
era is roughly a factor of two higher than in a single MOS 
camera at 2 keV, rising to a factor of four by 8 keV. The 
best-fitting spectral parameters obtained for separate fits 
to the PN and MOS spectra are summarised in Table |l|, 
where the errors are quoted at a 90% confidence level. The 
simple absorbed power-law model provides a very good fit 
to all three data sets. In addition the agreement between 
the spectral parameters derived independently from the 
PN, MOS 1 and MOS 2 datasets is excellent. 

The spectral variations revealed in Fig. || were inves- 
tigated by extracting spectra from the PN dataset from 



a set of annular regions centred on the point of peak sur- 
face brightness^. The resulting spectra were then fitted as 
before with a simple absorbed power-law model. Initially 
we allowed to vary as a free parameter between the 
different spectral datasets but, in the event, excellent fits 
were obtained with JVh fixed at the value obtained for 
the core region (Table Q). The derived photon spectral in- 
dex shows a steady increase with radius from a value of 
1.63 ±0.04 at the centre to 2.45 ±0.06 at the outer edge of 
the low surface brightness halo (Fig. ||). We note that the 
rate of spectral softening with radius will be slightly un- 
derestimated due to spreading caused by the instrument 
PSF (see Fig.||); however the impact on the derived spec- 
tral indices is likely to be small (AT < 0.1). Fig. ^ further 
illustrates the variation of the X-ray spectrum from the 
core through to the outer halo region. 

We next investigated whether there is any evidence for 
thermal emission from the halo component, focussing on 
the PN spectrum extracted from an outer (r = 96" — 144") 
annular region. As a first step we fitted a two compo- 
nent model, comprising a power-law and a (solar abun- 
dance, MEKAL) thermal component. The result was 
a modest improvement in the % 2 with respect to the 
power-law model (see Table |^) which the F-test shows to 
be marginally significant for two additional parameters. 
However, this putative (~ 1 keV) thermal component ac- 
counts for only ~ 1% of the observed flux. The lack of 
prominent line features in the X-ray spectrum of the outer 
halo is evident in Fig. ^| which shows the fitting residuals 
to a simple power-law spectral model. One possibility is 
that the line emission is suppressed due to the fact that 
shock-heated X-ray emitting gas in the halo of G21.5-0.9 
has not had time to reach a state of ionization equilibrium. 
When we fit a non-equilibrium ionization model (the NEI 
model in xspec) we obtain the spectral parameters listed 
in Table ||. The NEI model provides only a rather poor fit 
to the data when the column density is constrained to the 
value measured for the SNR core. A rather better result is 
obtained when Nn is allowed to vary freely (the last entry 
in Table |J) but nevertheless the fit is not as good as that 
obtained for the absorbed power-law model. 

We also examined the spectrum of the brightest part 
of the spur feature (located directly to the north of the 
core) but unfortunately the signal-to-noise ratio was too 
low to distinguish between thermal and power-law spectral 
forms. 

The total integrated X-ray flux of G21.5-0.9 within a 
radius of 144", corrected for the line-of-sight absorption, is 
6.1 x 10~ n erg cm~ 2 s _1 in the 2-10 keV band. Assuming 
a distance of 5 kpc based on HI absorption measurements 
(Davelaar, Smith & Becker 1986] ), this implies an X-ray 



luminosity of 1.8 x 10 35 erg s T . 

Finally we have briefly considered the spectral prop- 
erties of the south-west point source. We obtained an on- 
source spectrum from the PN dataset using an extraction 



2 In this analysis we excluded a 16" radius region centred on 
the south-western point source. 



Warwick et al.: XMM-Newton Observations of G21. 5-0.9 



5 



+ 



50 100 
Radius (arcsec) 



150 



Fig. 5. The observed variation in the photon spectral in- 
dex versus radius as measured in the PN data. 
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Fig. 6. Upper panel: The observed count rate spectra 
and best-fitting absorbed power-law model for the core 
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inner (48'' 



< 



< 96") halo and outer 



(96" < r < 144") halo regions (shown in blue, green and 
red respectively). Lower panel: The ratio of the predicted 
to the observed count rate spectra for the outer halo re- 
gion. 

Table 2. Best-fitting spectral parameters for the outer 
halo region of G21. 5-0.9 



Model 




r 6 


kT c 


n c t d 


X 2 /dof 
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356/230 
277/229 


a Column density in 


units of 10 22 


atom cm" 


-2 





Photon spectral index 
c Temperature in keV 
d Ionization timescale in units of 10 s cm" 
e Fixed parameter 



cell of 16" radius. A complementary background spec- 
trum was taken (on the same CCD) using a cell of the 
same dimension sited equidistant from the centre of G21.5- 
0.9. Spectral fitting indicated a relatively hard spectrum 
(r = 1.61 ± 0.34 or kT - 7 keV) and a column density of 
1.03 ± 0.22 x 10 22 cm -2 . As the latter is about a factor 
two lower than the derived TVh of G21.5-0.9, we can con- 
clude that this point source is a foreground object with 
respect to the SNR. If we estimated its distance as ~ 3 
kpc (i.e. roughly half the distance to the SNR in line with 
the ratio of Nh for the two sources) then the inferred 2-10 
keV luminosity is 4 x 10 32 erg s _1 . This point source is 
positionally coincident with the emission-line star SS397. 

5. X-ray Timing Analysis 

We have carried out a preliminary search for pulsed X-ray 
emission using EPIC data from the full observation (i.e. 
including periods of relatively high background rate) . The 
photon events in the 1-10 keV band from a region of 8" 
radius centred on the core of G21.5-0.9 were extracted for 
all three cameras. Next the event times were barycentric 
corrected using the standard SAS task "Barycen" and the 
latest available "Reconstructed Orbit" file. Power spectra 
were then calculated for the combined data from all three 
cameras and the PN dataset alone. No evidence was found 
for a significant periodic signal. In terms of the amplitude 
of an underlying sinusoidal signal, the detection threshold 
was ~ 5 x 10 -2 ct/s in the frequency range 5 x 10 -4 - 0.17 
Hz for the PN + MOS combination (the upper frequency 
limit being set by the frame time of the MOS CCDs). 
Corresponding values for the PN-only data were ~ 3 x 
10~ 2 ct/s in the frequency range 5 x 10 -4 - 2 Hz. In terms 
of the pulsed fraction (of the core emission) the limits are 
roughly 3.5% and 5.5% respectively. 

6. Discussion 

In a recent paper reporting Chandra observations, Slane 
et al. (|2000 ) consider the possibility that the low-surface 
brightness halo surrounding the central synchrotron neb- 
ula in G21.5-0.9 might be the shell formed from the ejecta 
and blastwave of the original supernova explosion. The 
present XMM-Newton observations help delineate the spa- 
tial extent and morphology of this component. Crucially 
the new spectral information demonstrates that the ex- 
tended halo has a spectral form devoid of any significant 
line features. If the halo emission is thermal then the lack 
of line emission implies that the plasma is far from ion- 
ization equilibrium. Our NEI modelling confirms a low 
ionization state with n c t w 3 x 10 8 cm" 3 s. The temper- 
ature derived for the continuum (bremsstrahlung) emis- 
sion is 4 — 5 keV which is rather hot, even for a very 
young SNR. Assuming the distan ce to G21.5-0.9 is 5 kpc 
(Davelaar, Smith & Becker 1986 ) and taking a maximum 



shell expansion velocity of 10,000 km s , we can esti- 
mate the (minimum) elapse time since the shock heating 
of the bulk of the X-ray emitting gas in the outer halo to 
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be ~ 100 yrs. Applying the above constraint on n e t then 
gives n c < 0.1 cm -3 . However, the observed X-ray lumi- 
nosity of the outer halo, if interpreted as bremsstrallung 
radiation, sets a requirement for a significantly higher elec- 
tron density, n c fh 0.5 cm -3 . There would appear therefore 
to be considerable difficulties in interpreting the extended 
X-ray halo in G21.5-0.9 as a thermal shell. 

There are known examples of SNR shells which exhibit 
non-thermal X-ray spectra due to local particle accelera- 



tion ( e.g. S N 1006, Koyama et al. |1995| ; G347.3-0.5, Slane 
et al. 1999). However, in the case of G21.5-0.9 the lack of 



any limb brightening, the smooth transition in spectral in- 
dex throughout the remnant and the remarkable circular 
symmetry, all suggest an interpretation of the outer halo 
as an extension of the central synchrotron nebula. We con- 
clude therefore, on the basis of the present observations, 
that G21.5-0.9 is a true Crab-like system rather than a 
composite (centre-filled plus thermal shell) object. 

The observed spectral softening with radius in G21.5- 
0.9 most likely reflects the impact of synchrotron radiation 
losses on very high energy electrons as they flow from the 
region of the termination shock around the pulsar to the 



be explained if our line of sight is reasonably aligned with 
the spin axis of the resident pulsar, so as to give a near 
face-on view of a putative inner torus. This hypothesis 
might also explain the non-detection of pulsed X-ray emis- 
sion from the core of G21.5-0.9, but is hard to reconcile 
with the axisymmetric structure seen in high frequency ra- 
dio maps and the near radial distribution of the magnetic 
field lines (Fiirst et al. 1988). In this setting the spur- like 
feature which emanates northwards from the core could a 
"plerionic" wisp associated with a complex inner magnetic 
field and/or jet structure. 
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average surface brightness of an underlying radio halo (be- 
tween r = 60" — 130") to be roughly 5 times this thresh- 
old. The inference is that very deep radio observations 
may well reveal the presence of an extended radio halo in 
G21.5-0.9. 

The differences in the spatial morphology of the Crab 
and G21.5-0.9 could be due to the viewing orientation. 
Specifically the near circular symmetry of G21.5-0.9 might 



